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THERMAL RADIATION TO A FLAT SURFACE 
ROTATING ABOUT AN ARBITRARY AXIS 
IN AN ELLIPTICAL EARTH ORBIT: APPLICATION 
TO SPIN-STABILIZED SATELLITES 


by 
Edward I. Powers 
Goddard Space Flight Center 


SUMMARY 


The derivation of total thermal radiation incident upon 
a flat plate rotating about an arbitrary axis is presented. 
The functional relationships between position in an elliptical 
earth orbit and direct solar radiation, earth-reflected solar 
radiation (albedo), and earth-emitted radiation (earthshine) 
are included. The equations have been programmed for the 
IBM 7090 digital computer, resulting in solutions which 
relate the incident radiation to spin axis orientation and 
orbital position. Several representative orbits for a typical 
geometrical configuration were analyzed and are presented 
as examples. 
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THERMAL RADIATION TO A FLAT SURFACE 
ROTATING ABOUT AN ARBITRARY AXIS 
IN AN ELLIPTICAL EARTH ORBIT: APPLICATION 
TO SPIN-STABILIZED SATELLITES 


by 
Edward I. Powers 
Goddard Space Flight Center 


INTRODUCTION 


The satisfactory operation of an artificial satellite depends upon maintaining the payload 
temperature within prescribed limits. For example, the standard batteries employed in present- 
day spacecraft generally restrict the temperature limits to 0° and 40°C. Often experiments lo- 
cated within the satellite structure further restrict this variation. 


The temperature level of a satellite may be determined by solving the instantaneous energy 
balance 


= dT 
P + Sa A. t diy tq. = gel, Tr qt ° 
where 
P = internal power dissipation, 
S = solar constant, 


a, = solar absorptance, 

ə instantaneous projected area for sunlight, 
qur, = reflected solar radiation (albedo), 
q., = earth-emitted radiation (earthshine), 


o = Stefan-Boltzmann constant, 


e = infrared emittance, 


A, = total surface area, 

T4 = mean fourth power, surface temperature, 
WC, = heat capacity of the satellite, 

dT ; : 

dt = time rate of change of satellite temperature. 


If the average orbital temperature is being computed, the last term is dropped. In this case all 
heat input terms represent integrated orbital values. 


It should be noted that the above equation, as-a representation of the entire satellite, is greatly 
oversimplified. The values for < and a, generally vary over the surface and great fluctuations in 
skin temperature may exist. In practice the thermal analysis consists of the development of a 
thermal model which represents a fine mesh of interconnected isothermal nodes. The appropriate 
relationship between nodes in regard to thermal conduction and radiation interchange must be 
established. 


The energy balance for each node of a thermal model may be written 


P: + ə A, + Fab | + Ges, t 2 Ens (E Te.) ə y N ed (EAT 
m m 


dT 


n 
= oeA, TE + (We) q ' 


n 


where 
C,„ = conductance between nodes n and m, 
E] = effective emissivity between nodes n and m, 
F,, = Shape factor-area product between nodes n and m. 


Since the major interest at present is to determine the satellite temperature level and not the 
gradients within, a discussion of the terms in the first equation is in order. 


For most passive controlled satellites P is relatively small compared with the total radiation 
input and does not have a significant effect on the satellite mean temperature. The magnitude of 
this effect, however, depends upon the « of the surface (e.g. a surface with a low absolute « may 
raise the internal temperature significantly because the skin has a limited capacity for 


re-radiation). 


The remaining three heat sources, direct solar heating, earth-reflected solar heating (albedo), 
and earth-emitted radiation (earthshine), represent the significant inputs to the satellite. It is 
apparent that an adequate thermal design is predicated upon a reasonably accurate knowledge of 
these thermal radiation inputs. The major source of energy—direct solar radiation—is fortunately 
the most accurately obtained. Since the sun's rays impinging upon the satellite are virtually 


parallel, the problem is simply one of determining the instantaneous orientation of each external 
face with respect to the solar vector. 


The calculation of earthshine is considerably less precise, requiring fundamental assumptions 
to reduce the complexity of computation. These include the assumptions that the earth is a dif- 
fusely emitting blackbody and that the surface temperature is uniform at 450”R. This permits 
direct calculation of energy input from all visible positions on the earth. Since all locations supply 
varying inputs, an integral equation must be solved to obtain the total incident energy. 


The albedo determination involves a similar integration. The earth in this case is assumed to 
be a diffusely reflecting sphere. In addition, the source intensity is a function of the satellite loca- 
tion relative to the sunlit portion of the earth. 


For a nonrotating satellite whose orientation is "fixed in space" the calculation of the thermal 
radiation fluxes is performed at any interval during the orbit. At each orbital positiona particular 
flat surface may or may not "see" the entire part of the earth's cap visible from the satellite. In 
the latter case integration for albedo and earthshine must exclude the shaded portion of the cap. 
Spin-stabilized satellites which rotate uniformly about the spin axis greatly complicate the analy- 
sis. The rotation around an arbitrary axis means in general that the heat fluxes vary, since the 
visible portion of the earth is a function of this rotation. 


The analysis presented herein is based upon the energy impinging upon a flat surface whose 
orientation is defined in vector notation (by the normal vector), and which is rotating about an 
arbitrary axis. This has a much broader application than may be realized at first. Although the 
obvious application is for spin-stabilized satellites, the results apply to any body of revolution. 
Here the interest lies in the variation of flux about the axis rather than an average value per spin. 
Thus, with the orientation of a single flat plate, the impinging fluxes on cylinders and cones are 
obtained. A sphere or a shape with a variable surface curvature along its axis requires several 
or many such plates, depending upon the precision required. It can be seen that the thermal 
radiation to the entire satellite surface may be found by simply considering a handful of appro- 
priately oriented flat plates. 


The purpose of this paper is to present, in general form, the derivation of the governing 
equations for the radiation energy sources as stated. The equations refer specifically to a flat 
surface rotating about an axis whose orientation is arbitrary. The dependence upon orbital posi- 
tion is included; the results of the numerical integration of these equations, encompassing a suit- 
able range of applicable orbits, is presented. 


COORDINATE SYSTEMS AND VECTOR REPRESENTATION 


The primary coordinate system is fixed with respect to the orbit (Figure 1). The XY plane 
lies in the plane of the orbit with the x axis coincident with the line connecting the center of the 
earth (the focus of the ellipse) and perigee. 


As the satellite traverses the orbit, the 
instantaneous location is defined by the angle 
a. "The altitude, therefore, may be determined 
at any instant by 


(P +R.) (1+ e) 
A(a) = l+ ecosa “R (1) 


where 


altitude, 


A(a) 


P 


altitude at perigee, 


LEGEND R 


O - EARTH CENTER x 
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radius of the earth, 


eccentricity of the orbit. 


Figure 1—Coordinate systems. The orientation of the unit vectors (Fig- 
ure 2) representing the normals to the flat 
plates, N, and the solar vector, S, are speci- 
fied by the following angles: 
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angle between the projection of N on 
the XY plane and the X axis, 
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angle between N and the Z axis, 
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angle between the projection of S on 
the XY plane and the X axis, 
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angle between S and the Z axis. 
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The vectorial representations in the fixed 
Figure 2—Orientation of normal and solar vectors. coordinate system are thus 


N = cosß siny i+ sinß siny j + cosy k (2) 
and 
S = cosy sind i + sinv sink j + cos k . (3) 


It is convenient to introduce an instantaneous coordinate system (X’ Y’ Z”) whose origin is fixed 
in the satellite (Figure 1). The X”Y” plane lies in the orbital plane with the x’ axis coincident with 
the line from the earth's center to the satellite. zZ’ and Z have the same orientation. 


The vectors N and S in the primed system are 


" 


N cos (B-a) siny i’ + sin (8-c) sinyj' + cosyk” , (4) 


S 


cos (v-a) sink i' + sin (v-a) sinAj' + cosAk' . (5) 


In a similar manner the expression for the spin axis vector A may be shovn to be 


A = cos (6-a) sing i' + sin (d-a) sinu j’ + cosy k' (6) 


Š = angle between the projection of A on the XY plane and the X axis, 


u. = angle between A and the Z axis. 
PLATE ROTATING ABOUT THE SPIN AXIS | 


For spin-stabilized satellites the normal 
vectors N, which represent the orientation of 
the exterior surfaces, rotate about the spin 
axis. In determining both instantaneous and 
average heat fluxes, the orientation of Nas a 
function of this rotation with respect to the 


primed coordinate system (X” Y' Z') must be 
known. This is accomplished by defining a 
third coordinate system X” Y” Z” (Figure 3). 
The X” axis is coincident with the spin axis A. 
Y" is defined sothat it liesin the plane formed 
by X” and an arbitrary fixed position of the 
normal vector N .* 


In Figure 3, the folloving terms may be 


Figure 3—Double-primed coordinate system. 


evaluated: 
cos T = cos (B-a) sin y cos (6-a) sinu 
+ sin (ß-a) siny sin (8- a) sinu + cosy COSH , (7) 
sin Tt = l-cos?7 , 
N = cost i” + sinr cosx j” + sinr sink k” (8) 


It can be seen that the angle between the spin axis andN, can be computed only in the 
first and second quadrants of the X”Y” plane. This has little effect since < varies from 0 to 27 
(a rotation). İH the vector N, lies in the third quadrant (of the X”Y” plane) the initial position for 
rotation remains in the second quadrant. | The objective now is to relate the double-primed unit 
vectors to the primed vectors in order to determine N in the primed system. 


*Since a complete rotation occurs, the orientation of N, merely indicates the starting point. 


The projections of j” and k in the primed coordinate system are found in the following 
manner: The value for the unit vector k” is 


; i' hi” + TK? 
k” = a A = Er J ; (9) 
| Ax Nİ VESES 
where 
g = sin (6-a) sing cosy - cosu sin (8-a) siny, 
h = cosu cos (ß-a) sin y - cos (6-a) sinu cosy , 
I = cos (8-a) sing sin (Bra) sin y 


- sin (8-a) sinu cos ($-a) siny 


Similarly 


| (Ax N.) x A Oi ae (10) 
where 
m = hcosu - l sin (98- a) sinu , 
n = f cos (5=a) sinu = gcosu , 


o = gsin (6-a) sinu - hcos (ë -o) sinu . 


N may now be determined in the primed coordinate system (X' Y' Z'), By rewriting the double- 
primed unit vectors, 


i” = pi' + qj' rk’ 
j” = si'+ ty? + uk” ; (11) 
k” = vi’ + wy’ + xk ' 
where 
p = cos (d-a) sinu , 
q = sin (ó—-o) sinp y 
r = COSH , 


A 1 nn, 
Ya? + h? +7? 


By substituting for i”, j”, and k” in Equation 8 


N ai’ + by" kek”, (12) 
where 
a = pcosT + ssinT cosk + vsinT sink , 
b = gcosr + tsin7 cosx + wsinT sink , 
c = rcosr + usin7 cos, + xsinT sink , 


where « is the angle of rotation. N may now be evaluated at any position during rotation through- 
out the orbit in the primed system (X' Y" Z”). 


DETERMINATION OF EARTHSHINE 


The general expression for the radiation intensity dq, from a diffusely emitting source dA, , 
impinging upon a unit area, is as follows: 


Icosw cos 7 dA, 


dq = en ee (13) 
where 
I = source intensity, 
w = angle between the line connecting the unit area with dA, and the normal to dA, (N,), 
n = angle between the line connecting dA, with the unit area and the normal to the unit area, 
D = distance between the unit area and dA, . 


VVith the application of this equation to 
a flat plate of unit area at an altitude Aa), 
the radiation intensity impinging upon the 
plate is 


oT. b iş cosw cos?) dA, 


GE TT > al, p? (14) 
where 
o = Stefan-Boltzmann constant, 
T: = mean black body temperature of 
the earth's surface, 
Figure 4—Earth's cap visible from satellite ($ = 0 o 1 = the heat flux leaving the source 
when coincident with Y”). dA 
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From the geometry of Figure 4 (vhich is similar to a sketch used by Katz”) and the vector notation 
presented, the terms in the equation may be defined more specifically: 


D = - Hi’ + R, sind cos¢ j' + R, sind sino k” , 


= [A(a) + R, (17 cos6)l i’ + R, sin cos¢ $” + R, sin sinó k” , 


P= [May + R (12 cos Rain"? | 
š R. sin 
go Qov tan” Ayr R, (17cos0) | ” 
D:N 
cosn = İDİ : 


= {- a fata) +R, (1- cos6)İ + b R sin6 cos¢ +c R, sind sing } : 


1 
[A(a) + R, (1- cos9)1? + (RÊ sin? 9) 


e 


*Katz, A. J., “Determination of Thermal Radiation Incident upon the Surfaces of an Earth Satellite in an Elliptical Orbit,” Grumman 
Aircraft Engineering Corp., Rept. XP 12.20, May 1960. 
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Equation 14 may now be written in the form 


R 


e 


ə də da 27 

= EI, Fe e a N R, sind 

qu = =: cos<® + tan Aa) + R, (1- cos 8) i 
0 


0 


ë {a [- A(a) + R (1- cos6)) + b R sind cos¢ + c R sind sing} 


1 3/2 
a 2 . 
Mm PR. (1-cos 9)) 2 + RÊ sin? >) By PLA (15) 


The limits of integration are such that the part of the earth cap visible from the instantaneous 
location of the satellite is included. For most practical cases, at least a portion of the cap is not 
visible from a plate because of its orientation. Attempts to define the appropriate integra- 
tion limits for such cases are extremely laborious. Since the equations cannot be solved 
without the aid of a high speed computer, an alternative approach is employed.* The numerical 
integration includes the entire visible part of the earth cap as stated, but the contributions of the 
elemental area that the plate does not see are deleted if the local value of cos” is negative. 


Equation 15 represents the flux for the instantaneous orientation of the plate. Interest also 
lies in the determination of the flux for a complete rotation of the plate about the satellite spin 
axis. The average value for q, is therefore 


27 
1 
I E = qe dx (16) 
0 


where « is the rotational angle. 


DETERMINATION OF ALBEDO 


The calculation for the albedo input to an orbiting plate is similar to that for the earthshine 
but involves additional basic assumptions. The reflectance of the earth depends on what is the 
visible surface—land, water, snow, cloud cover, etc. Until precise data is available which indi- 
cates the functional relationship between the albedo and the visible surface, an approximate mean 
value of 35 percent appears satisfactory. 


The reflected radiation is assumed to be diffuse. The local earth-reflected intensity varies 
with the cosine of the angle between the solar vector and the local normal. Since the input to the 
plate is a function of the source intensity of the visible elemental areas, the reflected radiation 
can be expressed by 


*Katz, A. J., op. cit. 


-1 R, | 
b cos İk + R. 27 * R. sin 
q = — cos 49 + tan Aa) + R. (1-cos8) j 


0 0 


° fa - A(a) + R. (1- cos6)) + b R sin cosd + c R. sin 0 sind) . 


3/2 
x... İs +R, ER 2 + RÊ J ull. (17) 
where 
S = solar constant (440 BTU/hr/1t?) 
alb = albedo factor (0.35) 
Y = angle between the solar vector S and the normal to dA,, N., 


Cos y may be obtained from: 


cosy = S-N 
cos (v-a) cos@ sind + sin (v-a) sinó cos¢ sink + cosà sin® sino. (18) 
The area of the earth that contributes to the albedo input is the surface common to the sunlit 


part of the earth and the part of the cap visible from the plate. The remaining area visible from 
the plate but receiving no sunlight contributes nothing to the heat input. This is accounted for in 


the numerical integration by deleting inputs when cos y is negative. 


The average albedo flux received by the plate per rotation about the spin axis is 


27 
1 
q, = ə | q, dk . (19) 


av. 
0 


CALCULATION OF DIRECT SOLAR FLUX 


The calculation of direct sunlight involves the determination of the instantaneous projected 


area of the plate with respect to the solar vector: 


= §(S-N) 


Isp 


(20) 


= s [a cos (v-a) sind + b sin (v-a) sind + c cosa] * 


where S is the solar constant. 


10 


Because of its orientation the plate may or may not be facing the sun during a rotation. Neg- 
ative values of S - N indicate that it is not facing the sun. 


The average flux per spin is 


to 
Asr = On asr dk . (21) 
0 


Because the orientation of the plate is fixed in space, the heat flux is constant throughout the sun- 
lit portion of the orbit. To determine whether or not the satellite is within the earth's shadow at 
any orbital position requires two simple checks.* If both of the following expressions are satisfied, 


the satellite received no direct input from the sun: 
cosQ < 0 , 


İR, + A(a)] + | sinQ| <R, , 


where 0 is the angle between the solar vector S and X’, 


cos Q = S:i', 


cos = cos (v-a) sind 


APPLICATION 


Hypothetical geometric configurations (Figure 5) have been chosen to illustrate the use 
of the equations. The external surfaces are represented by the indicated normal vectors. 


*Katz, A. J., op. cit. 


(a) (b) 


Figure 5—Geometrical configurations. 
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The three sketches represent three spin axis Table 1 
Mean Solar Heat Flux per Rotation” 


orientations. 
Or, Ve 
Circular orbits of 150, 550, 1050, 2000, and ö (degrees) | £ (degrees) (BTU/hr /ft2) 
5000 statute miles, in which the solar vector lies in 
the plane of the orbit (minimum sunlight), were con- 
sidered. Figures 6-35 indicate the orbital variation 
of earthshine and albedo.* These heat fluxes rep- 
resent mean integrated values per rotation at the 
instantaneous orbital position. The direct solar 
flux is constant in sunlight for a specific orientation. 
Appropriate values are presented in Table 1. 
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“The National Aeronautics and Space Administration... shall... 
provide for rhe widest practical appropriate dissemination of information 
concerning its activities and the results thereof . . . objectives being the 
expansion of human knowledge of phenomena in the atmosphere and space.” 


— NATIONAL ÁERONAUTICS AND Space Act or 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 


important, complete, and a lasting contribution to existing knowledge. 


TECHNICAL NOTES: Information less broad in scope but nevertheless 
of importance as a contribution to existing knowledge. 


TECHNICAL MEMORANDUMS: Information receiving limited distri- 


bution because of preliminary data, security classification, or other reasons. 


CONTRACTOR REPORTS: Technical information generated in con- 
nection with a NASA contract or grant and released under NASA auspices. 


TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 


TECHNICAL REPRINTS: Information derived from NASA activities 
and initially published in the form of journal articles or meeting papers. 


SPECIAL PUBLICATIONS: Information derived from or of value to 
NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Washington, D.C. 20546 


